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Q ■ In this talk I briefly comment on the conventional lattice gauge fixing adopting 

a critical, even though constructive, numerical point of view. 

1 Standard Landau 

On the lattice the usual procedure accepted to compute gauge dependent matrix 
elements is summarized in the following formula defining the expectation value 
of a gauge dependent operator O: 

°\ 

S i l P) = \\ dUO(U G )e~? s W : (1) 

d . Z J 

where S(U) is the Wilson lattice gauge invariant action, G is the gauge transfor- 



mation projecting the links in the Landau gauge 



> 

with the gauge rotation given by U^(x) = G(x)U^(x)G' (x + fi) and the gluon 
field defined on the lattice in the standard way: 



d^A^ = + periodic boundary condition (2) 



A^x) 



U,(x) - Ul(x) 
2iag 



(3) 

Traceless 



In the lattice gauge theories where the links belong to a compact group, the 
gauge fixing is necessary only in the case of the measure of a gauge dependent 
operator. An expression similar to ([]]) but without having gauge fixed the links, 
defines the expectation value of a gauge independent operator. Moreover, in a 
lattice simulation there is no need to compute the Faddeev-Popov determinant 
because the correct adjustment of the measure, necessary in the case of gauge 
fixing, is obtained rotating the links in the chosen gauge. Therefore the complex- 
ity of the ghost technique is replaced by the numerical evaluation of the gauge 
transformations. The price to pay is the large amount of computer time spent to 
obtain numerically the gauge transformations. From a numerical point of a view 
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the values of a gauge dependent operators strongly fluctuate around zero if the 
gauge has not been fixed. In the case of an imperfect or inadequate gauge fixing 
the measure of a gauge dependent operator is affected by additional fluctuations 
to be summed up to the intrinsic statistical noise. 

The necessary steps bringing to the computation of the integral (J^) can be 
described as follows: 

• A set of N thermalized configurations {[/} is generated with periodic bound- 
ary conditions according to the gauge invariant weight e~ Sw( - u ^; 

• For each {U} a numerical algorithm compute the gauge transformation G; 

• The expectation value of an operator is given by the mean value of the 
values taken by the operator on the gauge rotated configurations: 

{0) Latt = 1- J! 0{U G ) . (4) 

{con/} 

The gauge fixing algorithm is based on the minimization of a functional 
constructed in such a way that its extrema are the gauge fixing transformations 
corresponding to the gauge condition. The F standard form for the Landau gauge 
is: 

Fu[G] = -ReTr £ /7 M G( %) (5) 

fl,X 

and the transformations G for which |^ = rotate the links in the gauge d^A G = 
0. The algorithm sweeps all the lattice many times and it stops when a prefixed 
quality factor is reached. 

It is remarkable that the eq. (|5|) does not correspond to the natural discretiza- 
tion of the continuum functional 

F A [G] = -Tr J d 4 x (A G (x)A G (x)) = - (A G ,A G ) = -\\A G \\ 2 , (6) 

according to the lattice definition of the gluon field (H) but it differs from that by 
0(a) terms. The form in eq. (|^) is adopted not only for its simplicity but also 
because its minimization enforces the following discretized version of the gauge 
condition 

A G (x) = £ (A G (x) - A G (x - A)) = (7) 
fi=i 

where A^ must be related to the links by the standard definition (||). 

In order to study the approach to the minimum, two quantities are usually 
monitored. The first one is F[U G ] itself, which decreases monotonically and 
eventually reaches a plateau. The other one, denoted by 9, is defined as follows: 

Q G ^\- J29 G (x) = i Y. Tr [A G (x)(A G )t(x)],^ / d 4 x Tr(d,A G )\ (8) 



2 



where V is the lattice volume. 

The function 9 decreases (not strictly monotonically) approaching zero when 
reaches its minimum. The desired gauge fixing quality is determined stop- 
ping the computer code when 9 G has achieved a preassigned value close to zero. 
The choice of the gauge fixing quality is a delicate point in the case of a simulation 
with a large volume and a high number of thermalized configurations. Of course, 
the better is the gauge fixing quality, the more computer time is needed. More- 
over it is impossible to know before computing the gauge dependent correlation 
functions if the choice done is suitable. So that, the stopping 9 value is normally 
fixed on the basis of a practical compromise between the estimated computer time 
and the gauge fixing quality. Sometimes, in the case of calculations performed on 
computers with single precision floating point, the maximum gauge fixing quality 
is limited by a value of the order of the floating point zero: 6 ~ 10~ 7 , this value 
is usually enough to guarantee the stability of gauge dependent correlators. 

2 Lattice Gribov Copies 

On the lattice a conceptual and numerical difficulty connected with gauge fixing is 
the existence of many different minima of the functional F[U] . The different gauge 
transformations determined by different minima are not equivalent each other and 
can be labelled with the value of the functional F. Of course it is unthinkable 
to succeed in reaching numerically the absolute minimum. The search of the F 
minima is at least as difficult as to find the lowest state of energy of a spin glass 
system with hamiltonian F. So that the condition (|2[) does not fix the gauge in a 
complete way generating on the lattice a problem analogous to the Gribov copies 
in the continuum fl], [!|. However the analogy is only formal because it is not 
possible to establish a connection among continuum and lattice copies. Moreover 
it is also likely that (many) lattice Gribov copies are spurious solutions due to 
the discretization ||. 

Actually, it must be noted that the presence of Gribov copies in the lat- 
tice Landau gauge fixing is a very common phenomenon. For example, it has 
been shown that when the minimization algorithm includes the over-relaxation 
technique, varying the value of the over-relaxation parameter u different lattice 
Gribov copies |4]] are generated. Of course there is no correlation between the 
convergence rate and the value of F associated with the particular Gribov copy 
found. 

The numerical effects of lattice Gribov copies can be divided into two cate- 
gories: the distortion of a measurement and the lattice Gribov noise. The typical 
example of a distortion due to the existence of Gribov copies is the measure 
of the photon propagator in compact U(l) in the so called Coulomb phase. In 
this case the measure of the photon propagator as function of the momentum, 
performed using the gauge fixing in the standard way, was affected by a not 
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regular behaviour H. This problem was associated with the distortional effects 
due to the Gribov copies. In fact, after having chosen the gauge fixed configu- 
rations nearest to the minimum of the gauge functional, the photon propagator 
became a smooth momentum function. More recent studies || show the details 
of the Gribov copies dynamics and provide a practical procedure to eliminate 
their effects. It is interesting to note that in the case of the measure of the gluon 
propagator in SU(3) there is no signal in the literature about a similar problem 
(for a recent review see ref. 0). The numerical simulations are performed in the 
Landau gauge and the various authors claim that the effects of Gribov copies do 
not affect the measure. 

Anyway, in the normal case in which there is no distortion due to Gribov 
copies, there should be an increase of the numerical fluctuations due to the incom- 
plete gauge fixing associated with the copies. An attempt to study the properties 
of this noise has been done in ref. |J taking as an example the measurement 
of the lattice axial current Za- This quantity is particularly well suited to the 
study of the Gribov fluctuations, because it is a gauge independent quantity but 
it can be obtained from chiral Ward identities in two distinct ways: a gauge in- 
dependent one, which consists in taking the matrix elements between hadronic 
states without fixing the gauge in the simulation, and a gauge dependent one, 
which consists in taking the matrix elements between quark states in the Landau 
gauge. In the intermediate steps of the numerical computation, the second pro- 
cedure takes into account gauge dependent matrix elements potentially subjected 
to the Gribov noise. Hence, there is an explicitly gauge invariant estimate of Za 
which is free of Gribov noise and which can be directly compared to the gauge 
dependent, Gribov affected, estimate. 

The results of the analysis can be summarized in the following way: 

• there is a clear evidence of residual gauge freedom associated with lattice 
Gribov copies; 

• the lattice Gribov noise is not separable from the statistical uncertainty of 
the Monte Carlo method. 

The global effect is not dramatic because the Za value obtained with the 
gauge dependent methods (1.08(5)) is close to the gauge independent evaluation 
(1.06(6)) and the jacknife errors are comparable. 

3 Gluon Field Definition 

In order to impose the Landau gauge on the lattice it is necessary to define the 
gluon field in terms of the links. It is clear that the definition given in eq. (|3|) is 
far from unique and it cannot be preferred, from the first principles, to any other 
definition with analogous properties. Moreover, in the general field theoretical 
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Figure 1: Comparison of the matrix elements of (A^Aj}^) (crosses) and the rescaled 
(AiAi) ■ Cf(go) (open circles) as function of time for a set of 50 thermalized SU(3) 
configurations at (3 = 6.0 with a volume V ■ T = 8 3 • 16. The data have been slightly 
displaced in t for clarity, the errors are jacknife. 



framework any pair of operators differing from each other by irrelevant terms, 
i.e. formally equal up to terms of order a, will tend, to the same continuum 
operator, up to a constant. It is has been shown in ref. || that this feature is 
satisfied at the non-perturbative level in lattice QCD. In fact, different definitions 
of the gluon field, at the lattice level, give rise to Green's functions proportional 
to each other, thus guaranteeing the uniqueness of the renormalized continuum 
gluon field. The relation between the two definitions can be expressed up to 
0(a 2 ) terms in this way (|ITJ: 

A'^x)=C(g )A ll (x). (9) 

Therefore for a Green's functions insertions the following ratio is expected to be 
a constant 



(. . . A'^x) . . .) 



C(g ) ■ (10) 



(. . . Apix) . . .) 

This relation has been checked numerically on the lattice by measuring a set of 
Green functions related to the gluon propagator for SU (3) in the Landau gauge 
with periodic boundary conditions. In Fig. [I] the Green functions {A^^} and the 
rescaled one Cf(g )(AiAi) are shown, where 

(AiAi){t) = ^EE^(A(x,i)^(y,0)) (11) 

i x,y 
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Figure 2: Typical behaviour of 9 and 9 vs gauge fixing sweeps at (3 = 6.0 for a 
thermalized SU(3) configuration 8 3 • f6. 



and the operator (A^A'^t) is obtained replacing in the same form the alternative 
definition: 



traceless 



Aiag 



(12) 



The remarkable agreement between these two quantities confirms the proportion- 
ality shown in eq. 

From the numerical point of view, however, the various definitions are not 
interchangeable. In fact let me suppose to fix the gauge of a thermalized con- 
figuration stopping the gauge-fixing sweeps when 9 < 10~ 14 and then define 9 
as having the same functional form of 9, as in eq. (H), but with replaced by 
A^ given in eq. (|12|) The values of 9 and 9 during the minimization of F are 
reported in Fig. ^ for a typical thermalized configuration, as functions of the 
lattice sweeps of the numerical gauge-fixing algorithm. As clearly seen 9 does 
not follow the same decreasing behaviour as 9: after an initial decrease, 9 goes 
to a constant value, many orders of magnitude higher than the corresponding 
value of 9. This difference, already noted in ref. ||, between the behavior of 9 
and 9 could cast some doubts on the lattice gauge-fixing procedure and on the 
corresponding continuum limit of gauge dependent operators. 

This paradoxical situation is due to the fact that 9 is an operator and the 
lattice can attribute a value to it only after averaging it over the gauge fixed 
configurations of the thermalized set. Hence the comparison reported in Fig. [2] is 
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devoid of meaning because it is done comparing the values obtained by a single 
configuration. Moreover the behavior shown in Fig. can be readily understood 
in the following way. The operator 9, defined in eq. (^) (9 ), is computed in 
the lattice units taking the definition eq.(|) (eq.(|l^)) without the powers of a 
to the denominator. Then in the continuum variables 9 = ^1 d^x^d^A^x)) 2 
where V is the 4- volume in physical units (analogously for 9'). Hence, while 9 
vanishes configuration by configuration, as a consequence of the gauge fixing, 9 
is proportional to (d^A'^) 2 , which has the vacuum quantum numbers and mixes 
with the identity. The expectation value of (d^A'^) 2 , therefore, diverges as ^ so 
that 9' will stay finite, as a — > 0. 



4 Summary and Addendum 

Every step of the usual gauge fixing procedure is affected by subtleties. The 
Gribov copies can be moderately dangerous in a simulation but it is necessary 
to check their influence in any calculation. The definition of the gluon field in 
terms of the links is not a fixed prescription of the theory but it can be chosen, 
for example, in order to satisfy practical requests. 

It is also possible to take advantage from this freedom as it has been done in 
ref. [11, [12| in order to implement a procedure to fix a generic covariant gauge 



on the lattice. The great advantage of a covariant gauge is that varying the value 
of the gauge parameter it is possible to check numerically, in the calculation of 
gauge dependent Green's functions like for example the gluon propagator, the 
gauge dependence of the fitted parameters. 

After the end of this workshop a thorough study of the lattice covariant gauges 
and their applications has been completed |L3] . 
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